Theoretical rate of dissociative recombination of HCO+ and DCO+ ions 
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The article presents an improved theoretical description of dissociative recombination of HCO'*' 
and DCO'*' ions with a low energy electron. In the previous theoretical study (Phys. Rev. A 
74, 032707) on HCO^, the vibrational motion along the CO coordinate was neglected. Here, all 
vibrational degrees of freedom, including the CO stretch coordinate, are taken into account. The 
theoretical dissociative recombination cross-section obtained is similar to the previous theoretical 
result at low collision energies (<0.1 eV) but somewhat larger at higher (>0.1 eV) energies. There- 
fore, the present study suggests that motion along the CO coordinate does not play a significant 
role in the process at low collision energies. The theoretical cross-section is still approximately 2-3 
times lower than the data from a recent merged-beam experiment. 

PACS numbers: 34.80.Ht 34.80.Kw 34.80.Lx 
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The HCO+ ion has been known for more than a 
century and was the first ion discovered in interstellar 
space via microwave spectroscopy [l|. It was classified 
as unidentified element at the time of discovery. Later 
Klemperer \2\ suggested the linear HCO"*" ion as a can- 
didate, which was later confirmed by experiments (see, 
for example Ref. [3]). These last two decades, HCO+ 
and other small polyatomic ions have extensively been 
studied theoretically: their electronic structure, potential 
energy surfaces and equilibrium geometry have been sys- 
tematically investigated as, for example, in Refs. 0, Q. 
The spectroscopy of neutral HCO in the energy range rel- 
evant to these DR studies has been extensively mapped 
out by Grant and coworkers [6|. Dissociative recombi- 
nation (DR) of molecular ions like IICO+ plays an im- 
portant role in chemistry of interstellar clouds and there- 
fore allows astronomers to probe remotely various char- 
acteristics of these clouds. In space, the HCO+ ion can 
be formed by several possible associations such as II2 -I- 
C0+ or CH -h O and also H+ + CO [7, 8] and destroyed 
by DR. Different types of laboratory experiments have 
been performed in order to study DR of HCO+: after- 
low plas ma, merged-beam and storage ring experiments 
M, lly, [ill, \l2, [3 ■ From the laboratory experiments 
it is now known that DR in HCO^ proceeds mainly into 
the H+CO channel: HCO+ +0" ^ H + CO. On the 
other hand, at present, there is no consensus among dif- 
ferent experimental measurements of the actual DR rate 
coefficient: as they differ by up to a factor often [lll.ll2|. 

The theory of DR in diatomic ions has been reasonably 
well developed in recent decades. For triatomic ions, only 
recently has theory been able to provide meaningful re- 
sults for the simplest triatomic ion, HJ [14|, [islllql . The 
theoretical description of DR in triatomic molecular ions 
is a difficult problem in part because several different 
(electronic and vibrational) degrees of freedom have to 
be taken into account. Several approximations have also 
been made in a recent theoretical study of DR in HCO+ 
IITI. The obtained DR cross-section was about factor 2.5 



smaller than the lowest experimental cross-section [Tl| . 
One of the possible reasons why the theoretical cross- 
section was smaller in Ref. (171] than the experimental one 
is the approximation of the frozen CO coordinate fixed 
at its equilibrium value. Although the main dissociation 
pathway does not involve this coordinate, it was argued 
[17| that the CO vibration could increase the probability 
to capture the electron and increase the overall DR cross- 
section. In the present study, we improve the previously 
developed DR treatment in IICO+ and investigate the 
explicit role of the CO vibration. 



The theoretical treatment presented here resembles in 
many respects the approach applied previously to the 
HJ [lal and HCO"'" HTJ target ions. Below, we describe 
the new elements of the theoretical approach. We rep- 
resent the Hamiltonian of the ion-|-electron system as 
H = Hion + H^i , where Hion is the ionic Hamiltonian and 
Hf^i describes the electron-ion interaction. Consider first 
the ionic Hamiltonian Hion written in the center-of-mass 
reference frame. We use Jacobi coordinates to represent 
all vibrational degrees of freedom: Introducing G as the 
center of mass of C-0, the set of internuclear coordinates 
is represented by the quartet Q={RcOi Rgh, d, ip}. Here 
Rco and i?GH represent respectively the distances C-0 
and G-H, 9 is the bending angle between OC and GH, ip is 
the azimuthal orientation of the bending. Here, we con- 
sider i?GH as the adiabatic coordinate representing the 
dissociation path. In the previous study |l7| , the inter- 
nuclear distance Rco was fixed at its equilibrium value 
(-Rco=2.088 a.u) and the i?cH coordinate was treated as 
the dissociative coordinate. Note that even though we 
use an adiabatic representation, our inclusion of nonadi- 
abatic coupling effects means that we are not utilizing an 
adiabatic approximation since the adiabatic eigenstates 
form a complete basis set. 

The vibrational Hamiltonian in the Jacobi coordinates 



J^ion — 



2AiH-CO-RH-CO 



-co C^i^H-CO 



V{Rco,RGii,d), (1) 



where /xco and /in-co are respectively the reduced 
masses of the C-0 and H-CO pairs; L'^{d,(p) is the 
famihar operator of angular momentum corresponding 
to relative rotation of H and the CO axis. Represent- 
ing the ion by the above Hamiltonian depending on 
{Rco,Rgh,S,ip} only, we neglected by rotational mo- 
tion of the CO bond in space, but included relative ro- 
tation of H and CO. This approximation is justified by 
a large CO/H mass ratio. As a result of the approxi- 
mation, the projection m^ of the angular momentum L 
on the CO axis is conserved. We solve the Schrodinger 
equation with the Hamiltonian ([T|) keeping the Rgh coor- 
dinate fixed; this determines vibrational wave functions 
^m^.iiRGRi Rco,(^,'p) and corresponding adiabatic en- 
ergies Um^jiRGu) that depend parametrically on i?GH- 
Several of these curves are shown in Fig. [TJ The lowest 
adiabatic curves can approximately be characterized by 
quantum numbers {vi , Wj "^ , W3} of the four normal modes 
of the HCO+ ion. 
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FIG. 1: (Color online) Adiabatic curves Um^^i{RGu) versus 
the adiabatic coordinate Rgh- The curves are labeled with 
quantum numbers {vivl^'^ V3} of normal modes of HCO^. 
The «i quantum in our model corresponds approximately to 
motion along -Rgh and therefore is not defined for the adi- 
abatic curves because -Rgh is not quantized. The curves 
{1112 0} and {wi2 0} are almost degenerate. 
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TABLE I: Comparison of vibrational energies obtained in the 
adiabatic approximation with the exact calculation from Ref. 
[5|. The result of the previous study [131, where a different 
adiabatic approximation was used is also shown. In that study 
CO was not quantized and thus, 00° 1 was not calculated. The 
overall error is about 12%, which translates into about 25 % 
for vibrational wave functions. The energies are given in meV. 



table, since a similar (but not identical) adiabatic ap- 
proximation was used in that study. 

The structure of the electronic part H^i of the total 
Hamiltonian H is the same as in the previous study [l3| : 
R includes nsa, np'K~^, npa, npTr"*"^, and nda electronic 
states only. In the basis of the five electronic states, the 
Hamiltonian has the following block-diagonal form for 
each principal quantum number n: 
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(2) 



where Esa, Ep„, Epy^, and E^a are the energies of the 
corresponding electronic states at the linear ionic con- 
figuration; S and 7 are the real non-Born-Oppenhcimcr 
coupling parameters. The couplings S and 7 depend on 
-Rgh, RcOi and 9 and responsible for the Renner- Teller 
interaction. They are zero for linear geometry of the ion. 
The parameters in the above Hamiltonian are obtained 
from ab initio calculations of Ref. [i3| as discussed in 
Ref. [131 • In the present method, the electron- ion inter- 
action Hamiltonian -ffmt(Q) is now used to construct the 
5x5 reaction matrix Ki^i'{Q) written in the same basis 
of electronic states as Hint{Q)- 

Once the adiabatic states ^rn.^,i{RGii',Rco,S,f), en- 
ergies Urn ,;(i?GH), and the reaction matrix Ki^ii(Q) are 
obtained, we take .Rgh as the adiabatic coordinate and 
apply the quantum-defect approach that has been al- 
ready used in a number of DR studies of triatomic and 
diatomic ions [15|, llSl . We construct the reaction matrix 

^j,j'(^Gh) 



A reasonable measure of accuracy of the adiabatic ap- 
proximation is provided by comparing the energy split- 
ting between our adiabatic potentials Um^jiRcii) with 
exact calculations of the corresponding vibrational split- 
tings 5| (Table |T|) . The obtained values are practically 
the same as in the previous study ,17l] also shown in the 



K^{m^,l,i}.,{7n^,l.i}'iRGu) = (*m^,i|-^i,j'(Q)|$m;,i') , (3) 

where the integral is taken over the three internuclear 
coordinates, -RcOj f, and 9. The reaction matrix ICjj' 
thus obtained has many channels and parametrically de- 
pends on Rgh- For each -Rgh value, we then obtain a 
number of resonances with energies Ua{RGii) and widths 



Tal^GH)- The resonances correspond to the autoionizing 
electronic states of the neutral molecule at frozen Rgh- 
The fixed- i?GH width of the resonances is the reciprocal 
of the fixed- i?GH resonance autoionization lifetime, which 
of course is not an experimentally-observable resonance 
width since Rgh has not been quantized. 

The energies and widths of the resonances are then 
used to calculate the cross-section for electron capture 
by the ion. Depending on whether or not a particular 
neutral potential curve UaiRcn) is energetically open for 
direct dissociation, two different formulas are appropriate 
to use for the cross-section calculation. For the neutral 
states energetically open for direct dissociation, we have 



2^2 



U2 Z^ 
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IxURgh)? 



(4) 



However, the following equation should be used 

27r2 



^ = TtE l<X"^(^GH)|Vra(i?GH)|x+(i?GH))|'n3 



(5) 

for the UaiRan) curves that are energetically closed to 
direct dissociation [17[. In the above equations, ko is the 
asymptotic wave number of the incident electron, which 
depends on the initial state o of the target molecular 
ion; xt^R-GB) is the initial vibrational wave function of 
the ion; x^^^ {RqyC} is the quantized radial wave function 
of the UaiRcii) curve. The sum in the first equation 
includes all neutral states open for direct dissociation. 
The sum in the above equation is over all closed channels 
c that produce potential curves UaiRcn) closed to direct 



dissociation [17 1 



The total DR cross-section for HCO+ is mainly de- 
termined by the second sum because the electron is most 
likely captured into one of the lowest closed channels that 
cannot dissociate directly. As mentioned previously, the 
formulas above describe the cross-section for capture of 
the electron. It is equal to the DR cross-section only if 
the probability of subsequent autoionization is negligible 
compared to the dissociation probability, after the elec- 
tron has been captured by the ion. 

The projection M — m^ -(- A of the total angular mo- 
mentum on the CO molecular axis, where A is the pro- 
jection of the electronic angular momentum on the CO 
axis, is a conserved quantity in our model. Therefore, 
the resonances and the cross-section are calculated sepa- 
rately for each value of |A/|. Since A can only be or ±1 
in our model {a and tt states) and the initial vibrational 
state of the ion has m^ — 0, the possible values of \M\ 
are and 1. The total cross-section for electron capture 
by the ion in the ground vibrational level is given 17] by 



(cr*"*"') = (( 



,M=0 



) + 2(a^=i) 



Figures [21 and [3] summarize the results of the present 
calculation. Fig. [2] compares the present results with the 
experimental data from a merged-beam experiment |ll| 
and with the previous theoretical study [i3|. The theo- 
retical results are almost identical (about 10% different) 
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FIG. 2: (Color online) The figure shows the calculated DR 
cross-section for HCO^ (solid line) as a function of the inci- 
dent electron energy. The experimental [ll| (cross symbols) 
and previous theoretical [13] (dashed line) cross-sections are 
also shown for comparison. The theoretical curves include a 
convolution over the experimental electron energy distribu- 
tion according to the procedure described in Ref. [l9| with 



AEa 



A£|| = 3 meV. 



for electron energies below 0.1 eV. However they differ 
significantly at higher energies, where the present calcula- 
tion gives a higher cross-section. Both curves are smaller 
than the experimental data by a factor of 2-3. Therefore, 
the approximation of the frozen C-0 bond employed in 
Ref. [13] is apparently justified for low electron ener- 
gies but appears to deteriorate at higher energies. This 
result can be rationalized as follows. For small electron 
energies, the CO vibration plays a negligible role because 
only a few resonances are associated with excited CO vi- 
brational modes. In addition, normally, widths of these 
resonances are relatively small due to small relevant cou- 
pling in the corresponding reaction matrix elements: The 
largest coupling elements in the matrix are associated 
with the Renner- Teller coupling, which is active when 
rriip is changed. However, when the total energy of the 
system becomes close to (but below) the first CO-excited 
level {00°1} of the ion (0.3 eV above the ground vibra- 
tional level), the Rydberg series of resonances associated 
with the {00" 1} level becomes more dense and, more im- 
portantly, they become mixed with the Rydberg series 
of the resonances associated with {03^0}. The latter are 
coupled relatively strongly to the ground vibrational level 
{OO'^O} of the ion by the Renner- Teller coupling. 

Figure [3] shows the thermal rate coefficients obtained 
in the present study for HCO+ and DCO+ and compares 
them with available experimental data. Somewhat anal- 
ogous to Fig. [21 the theoretical DR rate coefficient for 
HCO''" is smaller than the rate coefficient obtained from 
the merged-beam experiment |11| by about a factor of 
3. The majority of the other experimental thermal rate 
coefficients shown in the figure are obtained in plasma 
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FIG. 3: (Color online) Theoretical (dashed lines) and exper- 
imental DR thermal rates for HCO+ and DCO+. The only 
available experimental data point for DCO^ is shown as a dia- 
mond symbol. The other symbols and the solid line represent 
data from experiments with HCO^. 



experiments. These rates are significantly higher than 
the merged-beam experimental data [ll[ . 

In the previous treatment ^^ , the DR rate coefhcient 
obtained for DCO"'" was approximately 30% smaller than 
for HCO"*-. In the present study, the DCO"*" coefheient 
is smaller than the one in HCO"*" by only 10%. The 
only experimental data available for 000+ is the rate 
coefficient 2.6 x IQ-'^cm^/s obtained at T=95 K ^. 
The same study provides the rate coefficient for HCO+, 
which is by 10% larger. 

Conclusion. The theoretical DR cross-section obtained 
in the previous study [17] was smaller by a factor 2- 
3 than the lowest measured experimental cross-section. 



However, the previous theory did not account for vibra- 
tion along the CO coordinate, and the main purpose of 
the present study was to assess the validity of the frozen 
CO approximation employed there. It suggests that the 
CO vibration does not play a significant role in the DR 
process at energies below 0.1 eV, but starts to be im- 
portant at higher energies, when the total energy of the 
ion-felectron system approaches that of the first excited 
CO vibrational mode. This study suggests also that re- 
duced dimensionality can be used in DR studies of small 
polyatomic ions as long as one includes (1) the dissocia- 
tive coordinate and (2) the vibrational coordinates re- 
sponsible for the highest probability of electron captur- 
ing. For HCO"*" the dominant dissociative coordinate is 
the CH bond (or Rgb), whereas the vibrational coordi- 
nates responsible for the electron capture are 6 and ip. 
Comparing with our previous study [i3| , the effect of the 
CO vibration seems to be larger for DCO"*", which can 
be explained by a larger D/CO mass ratio. 

In the present and previous theoretical studies, it was 
assumed that s- and p- wave-dominated eigenchannels are 
not mixed. This is justified to some extent by the fact 
that the ah initio energies used here account for the mix- 
ing, at least, at static geometries of the ion. However, in 
the dynamical framework of electron- ion collisions, the 
energy eigenstates obtained in the ah initio calculation 
could in principle be strongly mixed due to the perma- 
nent dipole moment (~3.5 debye) of HCO"*". The effect 
of the ionic dipole interaction with the electron should 
be addressed in future theoretical studies. 
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